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Abstract 

Various rescattering sources are analyzed in the context of the SU(3) flavor sym- 
metry. In particular, the possibility to account for intermediate charm at the hadronic 
level in B to PP is thoroughly investigated. Then, the rescattering sources are com- 
pared in light of recent B to two charmless pseudoscalar decay measurements, with 
emphasis on the size of strong phases and on patterns of direct CP-asymmetries. 
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1 Introduction 



Studies of CP violation in the B system offer a unique window into the intimate structure 
of the Standard Model, and possibly about New Physics. The dedicated experiments at 
Cornell, KEK and SLAC have begun to feed theorists with precise data to test the current 
understanding of CP violation, as well as alternate theories. 

Recently, a number of the rare hadronic decays B to two charmless pseudoscalars have 
been measured (see [1-6]). These modes are interesting to constrain the weak angle 7 
of the unitary triangle. However, for this program to be achieved, theoretical control is 
needed over the strong final state interactions (FSI). In the present work, this problem is 
addressed phenomenologically, at the hadronic level (for other approaches, see for example 
[2]). The basic tools are flavor symmetries like isospin or SU(3) (see for example [H], [0]) 
[TU] ) and Watson's theorem, or elasticity, for strong FSI (see ^0], f^)- 

Up to now, a serious shortcoming of elasticity-based approaches was their inability 
to account for intermediate charm. In short-distance analyses, quark diagrams involving 
charm quarks do have non-negligible imaginary parts ^Hl- This led to the common belief 
that rescattering cannot be elastic. One purpose of the present paper is to show that 
it is possible within the hadronic framework of Watson's theorem elasticity to account 
for intermediate charm. Some results in that direction were already presented earlier J2] . 
Here, we will revert to the language of the usual quark diagram topologies. The theoretical 
foundations of our parametrization, which are articulated around a factorization property 
for FSI, are reviewed in the first part of the paper. 

Another problem of elasticity-based approaches is that large strong phases are appar- 
ently required by current experimental data for B to two charmless pseudoscalars (see for 
example ^H). This is in disagreement with Regge computations |15j . with the expecta- 
tion for two particles flying apart with high momentum and with QCD-factorization 
predictions ^j. In the second part of the paper, we will show that when intermediate 
charm is accounted for, the strong phases tend to be much smaller. 

Finally, the last goal of the paper is to parametrize and quantify the corrections to 
the patterns for B ttK direct CP-asymmetries presented in (12j . These patterns were 
designed to discriminate among dominant rescattering sources, and we will here show that 
this capability is not altered. 

2 Theoretical Framework 

The starting point of our analysis is the "generalized Watson Theorem" (for details, see 
jlOj . jllj . jl2j). which permits the factorization of weak decay amplitudes into a factor 
invariant under CP and a factor that gets complex conjugated under CP 



W = ^/SWb 

CP (w) = Vsw^ 



(la) 
(lb) 
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where W is a vector containing the full decay amplitudes. Consequently, the real (up to 
CKM factors) Wb is identified as the bare decay amplitude (i.e. before FSI) denoted as 



W 



B ^ Xi 



(2) 



The unitary rescattering matrix v S contains all the CP-conserving strong phases and de- 
scribes final state interactions (FSI). In general, one defines the rescattering eigenchannels 



as the basis in which y 5" is diagonal 



Sdiag = Oc-S.O c 



e^^'^ca 



(4) 



with strong phases (or rescattering eigenphases) as diagonal elements. Elasticity 
(sometimes referred to as " quasi-elasticity" ) is defined as the conservation of the total 
decay probability, which follows from the unitarity of 



\W\ 



VsWb 



(5) 



Assuming that CP and CPT hold for the strong interactions generating v S implies that 
Oc is orthogonal (up to some phase conventions) since then is unitary and symmetric. 
The angles parametrizing Oc are then called mixing angles since they govern the FSI 
mixings among decay channels. 

In principle, y/S should be a n x n matrix with n the number of S-meson decay channels 
which can be connected by strong interactions. In practice, however, we expect that to a 
good approximation this matrix is bloc-diagonal. Unfortunately, it is not known at present 
how large the extent of each bloc should be to capture the essential rescattering physics. 
In the present work, we assume that two-body ^ n-body rescatterings are negligible (as 
could be justified from arguments ^H|)) or at least suppressed by large cancellations. 

The most restrictive (though non-trivial) approximation is that of SU(2)-elasticity: 
each bloc of corresponds to an isospin multiplet, which means that only rescatterings 
like vrvr ^ vrvr or KK ^ KK (see |12| ) are allowed. A more flexible approach is to 
consider SU(3)-elasticity, to open rescatterings between different iso-multiplets like vrvr ^ 
KK or ttK ^ r]sK. In an attempt to capture all the relevant physics, rescattering 
channels like DD ^ tttt or DgD ^ ttK should also be opened to account for intermediate 
charm. However, sizeable DD ^ PP rescatterings like those implied by SU(4)-elasticity 
would average decay amplitudes and produce Br (i? DD) ~ Br {B PP), in clear 
disagreement with experiment 



5r<=^P (B 



DsD) 



10^2 >> Sr^'^P {B K-k) ~ 10" 



(6) 
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Intermediate charm rescattering effects are thus quite small, and can be accounted for as 
distortions of SU(3)-elasticity. We call enlarged SU(3)-elasticity this parametrization of 
FSI, since the SU(3) flavor symmetry remains exact. Conversely, even if all the rescat- 
terings (both two-body and multi-body) are small, the DD ^ PP ones may still lead to 
sizeable effects because of Eq.©, so it is desirable to have an adequate formalism at hand. 

By definition, bare amplitudes W^, are real (except for CKM matrix elements) since 
they do not contain any FSI effects (see Eqs.(l)). To parametrize them, and to make 
contact with the short-distance weak decay of the b quark, it is convenient to introduce 
quark diagrams which account for the possible flavor flows (see, for example, 0,0) |Tn]). 
To first order in the weak interactions, and to all orders in the strong interactions (except 
FSI), the six topologies needed for B {PP}, P = 7r,K, % and B {DD}, D = D,Ds 
are shown in Fig.l. 



W 







Color-suppressed (C) 





Annihilation (A) Penguin (Pq) Penguin Annihilation (PAg) 

Figure 1: Quark diagrams for non-singlet final states, q = u,c,t. 

In principle, second-order weak interaction topologies are very small, except for those 
depicted in Fig. 2 involving an enhancement factor mf^p/M^. 




Electroweak Penguin [P^^) 




Color-suppressed 
Electroweak Penguin (^Pc^^ 



Figure 2: Second-order weak interaction quark diagrams. 

The SU(3) symmetry is implemented at the level of bare decay amplitudes by identifying 
quark diagrams that differ only by the interchange u ^ d ^ s. The badly broken SU(4) 
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symmetry is not invoked and, consequently, quark diagrams for B jDZ)} will be 
distinguished from B — > {PP} ones by a subscript D. 

Rescattering matrices have an interesting property of factorization. If each of the 
strong phases is expressed as the sum of two parts 

then, from the unitarity of (see Eq.(@}) 

Vsi5c,,...) = Vsi5^c,^...).VSi6'c„...) (8) 

Now, we can define effective " bare" amplitudes containing the part of the CP-conserving 
strong phases: 

W = \/5 (5c,, ...) . Wt = VS{6-c,,-) • (9) 

In just the same way Wi, is parametrized in terms of real bare quark diagrams T,C,..., 
the effective bare amplitudes W^^^ are decomposed into complex effective quark diagram 
amplitudes T^ff,C^ff,...{i.e. W^^^ = Wb{T ^ T^ff ,C ^ C^ff, ...)). It is important to 
realize that in the literature, quark diagrams are often introduced at the level of W^^^ . 
One should then be very careful in understanding which part of the FSI they contain, so 
that different approaches can be related. Finally, note that CP{W^^^) ^ [W^^^]* smce 
some CP-invariant strong phases are now included into them (compare with Eqs.(l)). 



2.1 SU(3)-elastic Rescattering 

SU(3)-elasticity is enforced by taking definite SU(3) states as rescattering eigenchannels. 
Taking as an example the set of decay amplitudes 



/ B+- 




\ 


B+ - 






B+ - 






\b^- 


DtD^ 


J 



(10) 



the rescattering matrix S is diagonal 

Sdiag = diag [e^''-\e'''^\e^''\e^''°) (11) 

in the basis defined by 

/ 127,3/2,1/2,1) \ 
127,1/2,1/2,1) 
18,1/2,1/2,1) 
V |8d, 1/2, 1/2,1) / 





3 

1 





30_ 
3_ 
10 



1 

V3 






1 / 



{K%+] 



(12) 
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Therefore, the SU(3)-elastic rescattering matrix is 

/ 



S 



SU{3) 



V 



9e'*27+e'*8 


v^(e»*27-e»''8) 


v^(e'*27-e*^8 


10 


10 

7e»^27+3e»'58 


3(6**8 -e'''27) 


10 


10 

3(e*''8-e*''27) 


^50 

2e'*27+3e»''8 




VEo 


5 












\ 



e< ) 



(13) 



See how SU(3)-elasticity fixes both the eigenphases Eq. l|llj) and the eigenchannels Eq.(|12|): 
the \D^D^^ and {-P-P} channels are not coupled and ^J^sui^) ^ characteristic bloc- 
diagonal form. Other sets of states are treated similarly (see !10|). 

Let us now concentrate on the PP ^ PP part of -^/^suis)- If each strong phase is 
expressed as 5j = 5- + 6f with i = 1,8,27 (one could think of the FSI short and long- 
distance parts), we can define effective, complex quark diagrams by absorbing the 6f part 
of the rescattering 

W^" = yJSsuis) iSh,SI,6i)Wb (14) 

Explicitly, effective quark diagrams can be expressed in terms of bare (real) ones as (omit- 
ting electroweak penguins for now) 



SU(3)-elasticity < 



x-ff 
ye// . 

peff . 



Pc,te'^^^ 



giA27-e'*8 



10 



(T + C) 
(T + C) 



X 
Y -- 



T,C 



E,A,Pu 



PAl^^ = PAub'^'^ 



PA 



eff 
c,t 



2Cl 



{T + C) + 
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(3r - C7 + 8^ + 8Pu) 



(15) 



The remaining 5' part of the rescattering is accounted for by acting with -^Ssuis) {^27^ <^8' "^i) 



on W^^^. 

In the following, we shall not use Eas-PSp. What we want to point out by writing them 
is that one should be careful when dealing with effective quark diagrams since obviously, 
the rough scalings (penguins are discussed in the next section) 



0{X) 



T 
C 

E,A 



(16) 



with A ~ 1/5 (21 may no longer be valid. For example, for the K^-k~^ mode, while A 

is helicity-suppressed and can be safely ignored, ^'^■^■^ receives unsuppressed contributions 
from T (provided (Jly — 5^ is non-negligible), as already emphasized in ^U]. Similarly, 
for — > vr+TT", E and PA can be neglected, but S*^'^-^ and PA^^^ may not. The 
Eqs. (|15j) illustrate in the specific SU(3)-elastic case a more general fact: if one starts with 
a parametrization in terms of effective quark diagrams without specifying anything about 
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FSI, all the effective quark diagrams are then arbitrary complex numbers and this costs 
much in terms of free parameters (except, of course, if one has a definite computation 
scheme for effective quark diagrams). 



2.2 Intermediate Charm Rescattering 



We now open the rescattering channel DD ^ PP |2Uj . For the set of states Ea. HlUj) . 
general rescattering eigenchannels are defined by mixing 8 with 8d (for a discussion on 
the formalism used here, see ^0], P^^: 

/Ci(27,3/2)\ /I 



C2 (27, 1/2) 
C3(8,l/2) 
V C4(8,l/2) J 



1 

COS X8 

\ -sinxs 
The rescattering matrix in the Ci basis is 














sin X8 
COSX8 / 



/ 127,3/2,1/2,1) \ 
127,1/2,1/2,1) 
18,1/2,1/2,1) 
V |8d, 1/2, 1/2,1) J 



diag 



(17) 



(18) 



since the mixing angle xs should be small (see Eq.®) 
then find 

/ 



In the physical state basis, we 



'5'sc/(3) +sin (2x8) 



— e 



2V10 



+ sin^ (X8) 
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-1 
V 






V 1 

V3 
-3 

-3^/2 








-V3 

-3^/2 

-6 









\ 





10 / 



1 \ 

-V3 
/ 



(19) 



with -^/Ssuis) given in Ea. ((T^ . In physical terms, one can understand xs as a DD- 
PP "coupling constant": the sin (2x8) term describes the DD pollution of PP states 
(which proceeds through DD ^ PP, DD ^ PP ^ DD ^ PP, etc) while the sin^ (xs) 
terms account for "inelastic" distortion of the rescattering among {PP} states due to the 
exchange of probability with the DD channel (PP ^ DD ^ PP, etc). The common factor 
g«<58 _ giSg g^g^g gg kiucmatical suppression since it decreases as the relative momentum 
between decay products increases 

(^8 = = in the limit Mb oo ^^). Finally, 
note well that SU(3) remains exact, and that the rescattering is elastic (conserved total 
probability) for the full set of states {PP, DD^ since ^J'S^ is unitary for all xs- 

Other set of states are treated similarly (the mixings among y = 0, Ts = states 
is described in Appendix A). All in all, six parameters are needed to describe all the 
rescattering mixings (one of the strong phase can be eliminated) 



PP : 



DD : 3 3 



27 



I X8 



1 

I Xi 
Id 



"8 '"1 



(20) 



6 



We expect that XSiXi << 1 since DD PP proceeds through the annihilation of the 
cc pair into hght quarks, and (^27, (^S)*^! < ^fi^f since there is less phase-space for DD 
than for PP. As said before, SU(4)-elasticity is not implemented: it would fix the Xi 
to some large values and introduce relations between 827, 6%, 5i and 5^,5i, leading to 
Br {B PP) K,Br[B^ DD) (see Eq.®). Let us stress again that SU(4) is used at no 
stage of the present analysis. 

The factorization property of ^/S can now be used to separate the DD ^ PP rescat- 
tering effects from the SU(3)-elastic part. For the set of states Ea. HlU|) . this is achieved 
by splitting the phases as (see Ec lfT^ ) 



(21) 



PP — PP 
DD ^ DD 



DD ^ PP 



The first factor does not depend on xs and is related to the SU(3)-elastic rescattering 
matrix Eo. 1)13(1 as S^g ((527, ^8) ^s) = \/^su(3) (^27,^8,^8), while the second factor con- 
tains all the effects of DD ^ PP mixing and becomes trivial if either 6g = 6s or X8 = 0- 
Proceeding similarly with the other multiplets, we now absorb the DD ^ PP part of the 
rescattering into effective quark diagrams 



On-shell cc FSI 



PAt" 



Pc-f3 



PAc + \f3 + a 



with 



1 



2\/3 



1 



,i(5f-<5i) 



XiTd, P 



XsTd 



(22) 



(23) 



To reach these expressions, we have retained only the dominant To contribution to B ^ 
DD decay amplitudes (see Fig.l) and the first order in Xi (so that PP — > DD — > ... — > PP 
effects are neglected, see Ea. lfT^ ). Importantly, all the other quark diagrams are unaffected 
by DD ^ PP. Since the Pc and PAc quark diagrams are precisely the ones involving the 
charm quark, we interpret Ea. l|22|) as a hadronic representation for on-shell intermediate 
cc quarks, as previously mentioned in jl2| (see Fig.3). 



B 





W 



B 



D 



D 



1::d p 



Figure 3: The standard quark- level penguin amplitude Pc and 
the hadronic- level penguin- like B ^DD^ PP contribution. 

In general, because of renormalization of the penguin loop, it is desirable to use the 
unitarity of the CKM matrix to write (Ag = V*ijVqj_(^s)) 

KPu + XcPc + XtPt = Xu (Pu - Pc) + Xt (Pt - Pc) = Xu (Pu-c) + Xt {Pt-c) (24) 
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At the i?-mass scale, we expect that Pu-c is small since mu^rnc « nib. Interestingly, 
when on-shell intermediate cc is treated at the hadronic level, the Pc short-distance ab- 
sorptive part has to be discarded and to a large extent the cancellation between P^ and 
Pc is preserved. Therefore, one can expect the rough scaling of the dominant penguin 
contributions [S] 



T 

Pt-c (25) 



pEW 



1 

0{X) 

and \Xu-c\ < \Xt-c\ with X = P,PA,P^^,P^^ to be valid. In the following, Xu-c wiU 
be neglected and the effects of intermediate charm will be accounted for by defining 

On-shell cc FSI | ^'^'-IfT ^'^'"^ " , , ^ (26) 

Note well that P^I^ and PA^[_^^ involve a combination of CP-conserving and CP-violating 
phases. 

Our method for implementing intermediate on-shell charm has semi- inclusive features. 
To see this, note first that all the developments above can be repeated replacing DD by 
D*D*,D**D**, ... Now, taking into account the effects on B —>^ PP of rescatterings from 
all the charmed meson states amounts to the definitions Eas. ()22() or Eas. H26() with a,f3 
replaced by 



a 



2^/3 " " V5 



P=JIfs, Fk= V (l-e'(4-S.))^^Ti (27) 



and this holds even if DD ^ D*D*, ... rescatterings are present thanks to the factorization 
property Eq.Q or Ea. (|21() . Finally, since in any case, the complex numbers 5 and (3 can 
be parametrized as in Eq. (|^ . this change is irrelevant (except, of course, that x and To 
can no longer be determined separately). In other words, as long as the Xi are small (i.e. 
when the last term of Eq. (|19j) can be dropped), Watson's theorem implies that no matter 
the precise physics in the charmed sector, its impact on i? — > PP will always amount 
to a redefinition of the penguin amplitudes Pc and PAc (or Pt-c and PAt-c). In the 
following, for simplicity, we assume that DD ^ PP dominates so as to give estimations 
of Xi- Finally, it is also straightforward to apply the formalism to charmonium modes Pr]c 
and/or rjoijc (see jU), simply by substituting Td by the Crj^ quark diagram contributing 
to Pr]c in P, and To by the C'^c contributing to tjqiJc in a. 

Let us conclude this section by a more technical comment. If we note 3, the bc.cd 
or bc.cs part of the weak Hamiltonian (see JHI)) Pc, PAc and the definitions Ea. ()22() are 
expressed in terms of SU(3) reduced matrix elements as 

' Pc = -(8|3e|3) f (8|3e|3)'^^ = (8|3,|3)+/3 

PA, = |(8|3c|3) + (l|3,|3) ^ I (l|3,|3)'=^^ = (l|3e|3) + a ^ ' 
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As expected, the a (/?) term involving xi (xs) contributes only to the reduced matrix 
element involving the 1 (8) final state, respectively. 



2.3 Final Parametrization 

All the pieces can now be put together to construct the final parametrizations of -B — > PP 
decay amplitudes assuming enlarged SU(3)-elasticity. From the previous sections: 

1- All PP ^ PP rescattering effects should be contained in \fS^ to preserve 
the scalings between quark diagram amplitudes, and will be treated assuming 
exact SU(3) (see Eqs.lO and Ens. flfill '). 

2- Intermediate charm can be treated at the hadronic level as DD ^ PP 
rescattering effects, and absorbed into Pt-c and PAt-c (see Ea.ipH))). 

The first point is important since it allows us to reduce the number of free parameters. 
Indeed, combining Eqs. (|l()j) and Eqs.(j2Sl) with CKM coefficients, it appears as sufficient 
to consider only four topologies (Ag = Vg*^V"gg') 

A5 = 0: \iT, XiPt^cXiC (29) 

AS" = 1 : AfPf_c, A^T, Af Pf_}^ 

since in any case, SU(3) breaking effects in the dominant T and Pt-c contributions should 
be greater than O (A^) ~ 4%. 

Once PA is neglected, it seems justified to set also (see Ea. ip^ ) 

2 \ - 5i = 5§ - 

/3 + a = 0^ ' V5 (30) 

since this combination describes rescatterings that proceeds through the vacuum (Eq.fPUJ 
may seem to imply a fine tuning between a and /3, but this is equivalent to the usual 
assumption that PAc is small, see Eas.ipS))). 
There remain eight free parameters 

QD amplitudes: T, C, Pt-c. PfJ^ 

SU(3) rescatterings: (527 — ^1-,^% — ^\ (31) 

On-sheU cc FSI: x^d, (5f - (5i 

The expressions for decay amplitudes are collected in Appendix B, in which the definition 
Ea. ()26() is used, b\ is set to zero and CKM coefficients are omitted. Also, the color-allowed 
electroweak penguin has to be introduced as C — > C + P^y[ [S]. Decay amplitudes for 
B DD, assuming that DD, PP ^ D*D* , ... are negligible, are also given. In that case, 
measurements of i? ^ DD branchings can serve to determine both To and the strong 
phase 5s — Si of PP ^ PP rescattering. This fact originates in the coherence requirement 
Eq. (|nn|) . which implies (provided X 7^ 0) 

Br (50 ^ D+D-) r (5°) _ 5 + 4cos {6^ - <5f ) _ 5 + 4cos (^g - ^i) 



Br{B+^D+DO)T{B+) 9 9 



(32) 
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3 Numerical Examples 



For our numerical examples, we consider only one source of rescattering: either SU(3)- 
elastic ones {PP ^ PP) or on-shell intermediate charm {DD ^ PP). The sets of free 
parameters in each case are 



Amplitudes: 
Rescattering: 



PP — PP 
brj - (5i, ^8 - ^1 



T T> C pEW 

Di 



DD ^ PP 



(33) 



with, in addition, the weak angle 7. To leading order, one can further neglect C and P^}^ ■ 
Though the formalism allows for general analyses, the precision of the experimental data 
is not yet sufficient for meaningful combined analyses. 

For each fit, we give the values of the parameters in the text, while the corresponding 
values for branchings (Table V), Acp (Table VI) and Sjrj (Table VII) are in Appendix C, 
along with the details of our fitting procedure. 

3.1 SU(3)-elasticity (no DD ^ PP) 

For the first example, we take the T and Pt^c amplitudes and — 5i, 827 — 5i strong 
phases (so we set X = 0)- For various (fixed) values of the weak angle 7 we find 



7 


60° 


80° 


100° 




18.3 


16.3 


17.7 


T 


0.65 


0.70 


0.76 


Pt-c 


0.11 


0.11 


0.10 


5s - 5i 


47° 


52° 


62° 


627 - 5i 


100° 


98° 


101° 



Table I: T,Pt_c, SU(3)-elasticity. 

Quark diagram amplitudes given in the table are dimensionless and produce branch- 
ings. 

Under this SU(3)-elastic parametrization, the B — > ttK branchings are dominated by 
the Pt-c penguin amplitude and verify 

Br {iT^K+) : Br (vr+A'^) : Br {k^K^) : Br {tt- K+) 1/2 : 1 : 1/2 : 1 (34) 

While the pattern of direct CP-asymmetries is (see |12j ) 



: A^oi^o : A^ 



1 3 

-K+ ^2: : : 1 

^ 2 2 



(35) 



no matter the strong phases. 

As already found in Jl], if the current CP-asymmetry measurement central values 
are to be trusted, they require large PP ^ PP phases. Then, in the B — > vrvr sector, 
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dominated by T, large rescatterings occur and the vr'^Tr" mode is fed from B — > {7r"'"7r }. 
The resulting pattern is close to 

Br (7r+7r°) : Br (vr+vr-) : Br (7r°7r°) 1 : 1 : 1/2 (36) 

Typically, Br {^B^ tt^it °) > 1.5xl0~6 under SU(3)-elasticity. For TT^TT time-dependent 
asymmetry, both A^+^- and 8^+^^- are found compatible with current experimental data 
(see appendix C): 




Y 



Figure 4: Fit result for A^+^- , S„+^- and a^fj under SU(3)-elasticity. 



A^+^- is not very sensitive to 7 and stays in the range 0.3 to 0.5. On the other hand, 
S^+^- varies much, and it appears that the present Belle measurement exclude 7 > 90° 
to 3a. However, the dependence of 3^+^^- on 7 is rather trivial. Defining the parameter 
Oieff as 

= yi^^42^sin (2ae//) (37) 

such that cteff = a = 7r — /? — 7if only T contributes to B^ vr+vr^ , the last graph shows 
that the deviation Oeff — ct stays small compared to 7. 

Finally, for B KK, no clear pattern emerges. However, one can note that because 
of large rescatterings from B — > {vr+vr"}, the B^ K^K~ branching is saturating its 
present experimental upper bound. Also, the B^ K^K^ and B^ — > K^K^ direct 
CP-asymmetries tend to be sizeable, between —1 and —0.5. 

Just for illustration, if we now introduce the subleading C and P^J^ topologies, the 
best-fit values are 



7 


60° 


80° 


100° 


A. mm 


10.1 


10.5 


15.4 


T 


0.60 


0.66 


0.67 


C 


0.21 


0.17 


0.17 


Pt-c 


0.11 


0.11 


0.10 


pEW 


0.011 


0.008 


0.005 


6s - 61 


53° 


57° 


32° 


627 - 61 


88° 


89° 


59° 



Table II: T,Pt-c,C,Pt^^ , SU(3)-elasticity. 

The most noticeable feature is a reduction of the strong phases, at the cost of a quite large 
C amplitude, and a significant decrease in Xmin- 
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To within 20%, the introduction of the additional quark diagram amphtudes does not 
modify the patterns described above. For example, the ratios of CP-asymmetries in the 
-kK sector are now 

8 13 3 , , 

A^OK+ ■ ^n+RO ■ ^nORO ■ A^-R+ = 2 + : -- + -K2/3 : -- + 3^0 : 1 (38) 

with, to leading order in |A^| T/ \ Xf \ Pt-c and P^}^ /Pt-c 

= COS 7 (cos (Jg - 527) (1 + ^) + a (1 - ^) ) (39) 



cos {68 - 827) - - 1 + 



Pt„c V 4 V C' + T 

For the range of parameter values in Table II, these corrections are of at most 15% (see 
Table VI). 

To close this section, it should be noted that the SU(3)-elastic constraint A^ox+ = 
2A^-j^+, though consistent, is not favored by the present measurements of A^-j^+ by 
Belle, BaBar and A^oj^+ by Belle: 

Belle BaBar 
A^-K+ -0.088 ± 0.035 ± 0.018 A -0.107 ± 0.041 ± O.OlSg] 
A^i>K+ 0.23±0.1lj:|j:[|^|21 -0.09 ±0.09 ±0.010 

Further, as shown in jl2j . this constraint is the same under SU(2)-elasticity. More precise 
measurements are necessary to draw any conclusion on the relevance of SU(N)-elastic 
rescatterings for B — > PP. 



3.2 On-shell cc rescattering (no PP ^ PP) 

For our second example, we take again the T and Pt-c amplitudes, but 5^ — 6^ as the only 
strong phase. In addition, we have the amplitude and the mixing angle x- As a first 
step, since it is always the combination (xTd) which appears in S — > PP decays, we can 
use the measurements .22 Br D' Df) = (0.8 ± 0.3) % and Br [B+ D^Dt) = 

(1.3 ± 0.4) % to fix Td = 2.43 which gives 

Q^th ^^0 ^ D-Dt) = Br*^' {B+ D^Df) = 0.98% (40) 
For various weak angle 7 we find 
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60° 


80° 


100° 


A. mm 


40 


32 


26 


T 


0.50 


0.57 


0.66 


Pt-c 


0.09 


0.10 


0.10 


X 


3.3° 


2.7° 


2.6° 


- Ss 


25° 


21° 


20° 



Table III: T,Pt-c, on-sheU cc. 
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Note that the fitting procedure is dehcate because, to a large extend, Pt-c and x^d 
compete against each other (see Eq.(j2ni)) and the Xmin function is rather flat. The above 
results are those for which both x and 6^ — are simultaneously small. 

The B ttK branchings are now dominated by both Pt-c and B — > ttK 
contributions, and since these effects behave like a "hadronic" penguin, see Eq. (|2(ij) . the 
pattern Ea. ()34() is not altered. Concerning CP-asymmetries, we find the pattern 

Pt Pt 

A^OK+ ■■ A^+KO : A^OKO : A^-j^+ ^ 1 : — : — : 1 (41) 

Interestingly, some ratios of CP-asymmetries directly give the value of Pt~c/T. These 
asymmetries are generated through the interference of T,Pt_c with B ^DDg^ irK, 
and require only one strong phase 5g—6s of less than about 25° (generated by the mismatch 
between PP ^ PP and DD ^ DD, see En.(l2ni)V 

In the B ^ TTTT sector, the fit is not very good (the large Xmin essentially generated 
by the and TT^TT^ modes). Indeed, the B — > ji^I?} vrvr effects are Cabibbo- 

suppressed, and since there is no averaging from PP ^ PP strong phases, the decay 
branchings keep their bare scalings in terms of dominant T (see appendix B) 

Br (7r+7r°) : Br (vr+vr") : Br (vrV^) w 1/2 : 1 : (42) 

which is rather far from current measurements (but not yet ruled out, see table V, appendix 
C). Typically, Br (i? — > vr^^vr'') < 10~^ without PP ^ PP strong phases (and assuming 
C/T ~ O (A) or smaller). Finally, A^+^- is found a bit smaller (in the range 0.2 to 0.4), 
while S^+^- behave roughly as in the SU(3)-elastic case (see Fig. 5). 




Y 



Figure 5: Fit result for A^+^~, S^+jr- and a^/f under on-shell cc rescattering. 



Note, however, that Oe// 7^ a is now due to both B \^DD^ — > tttt and Pt~c, and 
therefore that no Pt-c penguin does not imply cteff = ot (this is obvious from Eq. (|26j) ). 

For B KK, the branchings are similar than in the SU(3)-elastic case except for 
B~^ K^K~~ which is exactly zero (no rescattering from B — > {7r"'"7r~} and no B ^ 
{DD} K+K~ because of Eq.dSOl)). The B+ K+K^ and B^ K^K^ direct CP- 
asymmetries tend to be much smaller (between and —0.3), and in addition verifies 

^K+Ro '■ Aj^oj^o : A^o^^o = 1:1:1 (43) 

If we now include the subleading C and P^}^ , the best fit procedure produces a large 
C amplitude to match the observed pattern of -B — > vrvr branchings (we have restricted 
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C/T to be at most 1/2) 
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0.004 
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The pattern of B 
amplitude 



Table iV: l , J-'t-c,'-^', J-"t^c , on-siieii cc. 
ttK direct CP-asymmetries is much affected by such a large C 



A. 



+ K0 



A. 



A^ 



-K+ 



-, c 

IH 

T 



CPt- 

rp2 



Pt-c + C CPt^ 



T 



rp2 



1 (44) 



while corrections due to P^}^ are subleading. Notice, however, that for the range of values 
in Table IV, this pattern is still different from the SU(3)-elastic one Ea.(|5^. Finally, the 
pattern Ea. (|43|) survive only partially: while Aj^+f^o and Aj^oko remain exactly equal (no 
C or P^l^ contribution), ^7^0,^0 is dominated by C and varies much. 

To close this section, let us make a few comments. First, note that phenomenologically 
speaking, Pt-c/T = is not yet ruled out by experiment once the penguin-like rescatter- 
ing effects B — > {DZ)} — > PP are allowed. Indeed, the B Kit branchings could be 
dominated by the B {D^I)} — > ttK contributions, while both B ttK and B vr+vr^ 
direct CP-asymmetries could be generated by interference of B —> {DZ)} — > PP contri- 
butions with T. In other words, assuming that DD ^ PP is the dominant rescattering 

process, the current data requires P^ll /T ^ with P^ll given in Ea. H26() . but do not say 
anything yet about the relative strength of the quark-level penguin Pt-c vs. its hadronic- 
level counterpart B {Z)Z)} PP. Now theoretically, from short-distance analyses 
there is good reasons to expect Pt-c/T ^ 0, and we hope that the phenomenological 
patterns of CP-asymmetries described above will help in extracting it from the data. 

Second, under exact flavor SU(3), it appears that consistency with B — > -kK direct 
CP-asymmetry measurements restricts j4^^ to be smaller than ~ 0.5. This statement is 
independent of the precise physics in the charmed meson sector. Except for the values 
quoted for the results of the fits are equally valid when the intermediate charmed meson 
state is DD, D*D* , ... or a charmonium state like Pijc (see \21\). or also combinations of 
them. Further, the rescatterings between these states is irrelevant for B PP. Note also 
that abandoning the coherence requirement Eg. 1)301) does not lead to larger Aj^t^, even if 
PA^J^ = PAc + 2/3P + a ^ 2/2,(3 + a (see Eq.(I22I)) contributes to B^ tt+tt- and not to 
B — > ttK, because PA^J^ also contributes to B^ — > K^K~ and therefore cannot be large. 

Finally, it should be noted that many relations between B^,B^ direct CP-asymmetry 
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and Bs ones exist 



A 
A 



Bs^-K+K ^B0-»7r+7r 



^Bs^K+K- 



'- ^BO^n~R+ 



All these relations are unaffected by C or P^r}^ contributions and specific to the dominant 
on-shell cc rescattering source (no PP ^ PP rescatterings) . 

4 Conclusions 

The B to two charmless pseudoscalar decay processes are analyzed assuming exact SU(3) 
for the amplitudes and rescatterings. The basic tool is Watson's theorem for hadronic 
final state interactions, which allows for a factorization of rescattering effects from bare 
decay amplitudes. In addition, the various FSI processes can also be factorized and treated 
separately. In the present work, we assume SU(3)-elasticity for the PP ^ PP part of 
the rescatterings, while DD ^ PP elastic rescatterings are shown to amount to the 
redefinition of the penguin amplitude as 



where is the tree diagram contributing to — > DD and x the small mixing parameter 
governing the DD pollution of i? ^ PP (this parameter is not fixed by SU(3), and must 
be small because of Eq.Q). Importantly, the redefinition Ea. ()45() also apply when the 
intermediate charmed meson state is D*D* , D**]j** ^ Prjc, or any combination of them, 
and therefore Ea. (|45jl can be seen as a hadronic representation for on-shell intermediate 
cc. 

Also, it should be clear that once the various rescattering effects can be treated sepa- 
rately, different theoretical tools can be combined. This is well illustrated by Ea. ()45() : if one 
has a definite computation scheme for P^I^ accounting for (elastic or inelastic) intermedi- 
ate on-shell cc, but suspects non-negligible elastic long-distance PP ^ PP rescatterings, 
the present formalism permits the combination of the two. In other words, it is perfectly 
consistent to account for SU(N)-elastic long distance rescatterings by acting with the cor- 
responding y/Ssu{N) on effective quark diagram amplitudes already containing various 
rescattering effects, like those produced by QCD-based approaches. 

In our phenomenological approach, many amplitudes can be safely ignored and we are 
left with the set of free parameters Eg. pif) . with, in addition, the weak angle 7. To further 
restrict the number of free parameters, one rescattering source is assumed dominant: either 
SU(3)-elasticity (PP ^ PP) or on-shell cc {DD ^ PP). The most prominent features of 
each fit are: 

SU(3)-elasticity: The B ttK direct CP-asymmetries tend to require large 
strong phases. This is nice for B vrvr (see Ea. ()36() ) but may be problematic 
for B — > KK (those modes are predicted close to their present experimental 



X,P^l = A,P,„. - A.y| (1 - e^(^f -^«)) xTn 



(45) 
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upper bounds). Also, it should be noted that the SU(3)-elastic pattern of 
B — > ttK direct CP-asymmetries Ea. ()35() is rather insensitive to subleading 
quark diagram amplitudes (C and P^}^), and is not in good agreement with 
present measurements (though the discrepancy is not yet significant). 

On-shell cc: A single strong phase difference {6g — 6s) of around 25° is suf- 
ficient to produce the required asymmetries in the ttK and tttt sectors. This 
strong phase difference is generated by the mismatch between PP ^ PP and 
DD ^ DD elastic rescattering channels. Such a small strong phase is nice 
from the point of view of Regge computations and from the naive expec- 
tation for two particle flying apart with large momentum but also with 
the SU(2)-elastic isospin analysis of -B — > Dtt P3I . Compared with the SU(3)- 
elastic case, the fit is not very good in the B tttt sector (no rescattering 
from B vr^^TT" leads to the approximate pattern Ea. H42() l. but better for 
B — > ttK direct CP-asymmetries (the pattern Eq. (|4H) is favored by current 
measurements) . 

Combining both sources of rescattering would certainly improve the quality of the fit 
as more free parameters are available. For example, the difficulty of the on-shell cc fit 
in the B tttt sector is reduced by including a small amount of PP ^ PP effects (like 
^27 — Ss ^ 20° ) . More precise measurements are necessary to pursue the analysis in that 
direction. 
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A Rescattering Matrix for = Y = states 



This appendix presents the construction of the enlarged SU(3) rescattering matrix for the 
Ta = y = states. We start with 



127,2,0,0) \ 

127,1,0,0) 
127,0,0,0) 

|8, 1,0,0) 

|8, 0,0,0) 
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where (entries are written according to x = sign (x) ^/\x 


3/5 


2/5 





o 



pp 

SU{3) 



^SU{3) 



/ 










-1/5 


1/5 




-3/20 


-3/20 27/40 




3/10 


-3/10 




-1/10 


-1/10 -1/5 


V 


-1/4 


-1/4 -1/8 




1/2 


1/2 \ 


( 


-1/6 


1/6 2/3 




-1/3 


1/3 -1/3 / 



-1/3 


1/60 


2/5 
-1/4 



2/3 \ 


1/120 

1/5 
-1/8 J 



To introduce rescattering between DD and PP, it suffices to mix |8, 1, 0, 0) with |8d, 1, 0, 0), 
|8, 0,0,0) with |8d, 0,0,0) and |1, 0,0,0) with |1d, 0,0,0) by defining 
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Finally, the rescattering matrix is found as 
when the Xi are small, 



•^x = Og^,^yO^^.^/Sdi^.O^.Osu{3), with. 



mg 



diag ^ 



■2iS27 „2iS27 „2iS27 J2iSs 2iSs 2i5i 2iS^ 2i5f 2iS- 
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B Decay Amplitude Expressions 



The B PP decay amplitudes, under the approximations discussed in the text, are 

=f^{T + C)+e^'^Pt[i 



A{B+ ^ K+K^) = i^l^flL (T + C) + e'^^ P^^^ 



A{B^ ^ K+K-) = ^iizfzmi^zfliwr - ^^^-" '^-y — + '> p^li 

/lrRO^^+^-^_ 5+8e'^8+7e'^27 ^ 5(l-e'^8 )+2l(e'^8 ) 2+e'^8 pe// 

irBO^^O^O^_ 5(l-e'^8)+l3(e'^8-e'^27) 5+16e%+39e!^ ^ , pe// 
^ TT TT j - ^ 1 ^ G + -^Pt-c 

A{B+^K+n'^)= e-^^03L^T + C) + ^^P:Ii 
AiB+ ^ K°7r+) = ^^1^^ (T + C) + e'^^Pt'I[ 

A{B^ - K+TT-) = 3-"^+^-"" r + + e^^spe// 

A{B^^^-^n,)= ^^^^|^(r + c) + ^p,!f/ 

The ampUtudes for the charge-conjugate modes are identical, except for the understood 
CKM coefficients. Decay amplitudes for Bg mode are: 

A {Bs - K-^+) = ^A^T - ?(f!!!zf!!!!) c + e^^^PtU 

A(R. ^ K^ir^) - 3(e'^8-e'^27) 2e'^8 +3e'^27 ^ e'^8 pe// 



20 ^ 60 

+(e'^8-e'^27) _ 5(l-e'^8)-3(e'^8-e'^27) 2+e'^8 pe// 

20 60 3 t-c 

A{Bs ^ TT+Vr-) = 5(l-e'%) + (e^^8-e^^27) ^ _ 5(l-e'^8 )-3(e'^8 -e'^27 ) ^ ^ 2(l-e'^8) ^e, 



20V2 60a/2 3t/2 



A(B,^V^m)- 5+8e-^8+27e-^27 ^ 2(l+e'^8) 

^ ^ - 20V2 60V2 ^ ^ 3V2 < 



t-c 
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Finally, assuming negligible PP, DD ^ D*D*, D**D** , the B DD decay amplitudes 
are 



A{B+ - 




= a{bo- 
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- = 










a{bO- 


D+D~) = 


le^'s (2 + e- 




a{bo- 


> D^D^) = 


-A {B^ - 


- DtD~) = 


^e'^s (1 - e- 


-iSs^ Td 






A{B's- 


- DtD~) = 


le^'s (2 + e- 


A{Bs- 




-A{Bs- 


-^D+D-) = 


le'^s (1 - e- 



C Fit Results for Branchings and Asymmetries 

The experimental data on B to two pseudoscalar decays we shall use are summarized in 
table V and VI, where direct CP-asymmetries are defined according to the sign convention 

_ T{B^f)-T{B^f) 

- T{B^f)+TiB^f) ^''^ 

Average is made over CLEO Belle j2] and BaBar |3| measurements (see also assum- 
ing no correlation. The implications of SU(3) symmetry bear on the decay amplitudes, so 
predictions for branchings must be corrected to account for lifetime differences and avail- 
able phase-space. From |^, the lifetime correction factor is Ftot iB~^) /Tfot = 0.92, 
while that of phase space are of at most a few percents, and are neglected. The branch- 
ings under brackets in table V are branching fractions for the physical ij state and not for 
778 1 so that they are only indicative. Indeed, the small admixture of singlet tjq into the 
physical rj state could lead to large effects. The last observables of interest to us are the 
time-dependent asymmetry parameters 

|A..-|^-1 2ImA.r A(B^ff) 

The experimental situation and averages are (with inflated errors, see |22j): 

Belle BaBar jS] Average 

A^+^- 0.77 ±0.27 ±0.08 0.19 ± 0.19 ± 0.05 0.38 ± 0.27 
S^+^- -1.23 ±0.41 ±0.08 -0.40 ±0.22 ±0.03 -0.58 ± 0.34 

We have used only A^+^~ = —C^+^- as input. 

The function we minimize to find the best-fit values is as usual 

f \ '^Br J f \ ^AcP J 

where the sum runs over measured decay branchings and asymmetries (12 inputs). Up- 
per bounds are implemented using the arc-tangent representation of the step function (3 
inputs). 
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Finally, the Wolfenstein parametrization j24j of the CKM matrix elements is used 

VcKM = -A 1 - ¥ - (f + 1) (49) 

\ AA3 (1 - p - ir^) -AX^ + AX^ {\ - p - ir]) 1 - ^ / 

and we take the central values of [22] 

A = 0.2196 ± 0.0023, A = 0.854 ± 0.045, y/p^T^ = 0.43 ± 0.07 (50) 

and keep only the weak angle 7 = arctanr//p as a free parameter (so unitarity of Vckm 
is implied). For the parameter we use the average sin 2/3 = 0.736 it 0.049 extracted 
from charmonium modes (25) . 

The result for branchings and asymmetries for the fits are collected below. Column 
labels refer to the best-fit parameter tables given in the text. 
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[4.0±0.9] 

<1.3 
5.3±0.8 


1.59 1.38 
1.09 1.28 
3.48 4.42 


2.41 1.72 
1.01 1.28 
5.43 5.67 


1.42 1.17 
0.88 1.21 
2.07 3.10 


2.31 1.93 
0.89 1.18 
4.34 5.86 


B" 


K+K- 


<0.6 
<1.6 

[<2.9] 
4.6±0.4 
1.92±0.44 


0.60 0.60 
1.02 1.31 
0.17 0.25 
0.26 0.38 
4.82 4.79 
2.12 1.70 


0.60 0.15 
0.93 1.14 
0.12 0.09 
0.26 0.38 
4.58 4.53 
1.77 1.57 



0.81 1.11 
0.05 0.06 
0.27 0.37 
5.03 4.98 
0.41 0.56 



0.82 1.08 
0.11 0.09 
0.27 0.36 
4.80 4.47 
0.68 1.43 


AS = 1 


Exp. 


60° 90° 


60° 90° 


60° 90° 


60° 90° 


B" 


K+TT- 


[2.4±0.9] 
11.2±1.4 
18.2±0.8 


3.33 3.17 
9.98 9.51 
18.6 19.2 


2.94 2.82 
8.81 8.45 
18.3 19.0 


3.37 3.17 
10.1 9.51 
18.6 19.2 


2.92 2.86 
8.75 8.57 
18.5 19.1 


B+ 


K+V8 

K+TT° 
K°7T+ 


[3.1±0.7] 
12.8±1.1 
20.6±1.4 


4.02 3.52 
10.3 10.4 
22.0 20.6 


3.68 3.20 
11.7 11.7 
22.3 20.7 


4.02 3.51 
10.1 10.4 
22.0 20.6 


3.50 3.20 
11.6 11.6 
21.8 20.5 


Bs 


K+K- 
VsVs 




15.5 14.7 

16.6 14.5 
14.1 12.5 
0.03 0.02 
5.19 6.95 
2.59 3.48 


14.2 16.9 
15.9 17.0 
12.9 14.2 
0.05 0.04 
6.55 1.95 
3.28 0.97 


17.6 18.2 
19.2 18.0 
17.0 16.0 





17.6 18.1 
19.0 17.9 
15.8 15.2 
0.05 0.03 





Table V: Fit results for charge-average branching fractions. 
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Direct CP-asymmetry 


SU(3)-elastic 


On-shell cc 




AcP (%) 


I 


II 


III 


IV 


AS' = 


Exp . 


U I J 


uu 


ou 


(1(1 


OU 


(1( V "" 


in) 


uu 






— 


-78 


-99 


-88 


-70 


-32 


-19 


-72 


-36 


Bs 




— 


-78 


-99 


-88 


-70 


-32 


-19 


-72 


-36 




K-TT+ 


— 


23 


25 


19 


26 


23 


22 


26 


31 






— 


49 


60 


28 


41 


23 


27 


27 


35 


B+ 




— 


-71 


-64 


-67 


-56 


-32 


-19 


-33 


-24 






-7±15 




























K+K- 


- 


14 


16 


11 


18 


- 


- 


- 


- 






— 


-88 


-95 


-87 


-58 


-32 


-19 


-33 


-24 




VsVs 


— 


-96 


-98 


-16 


39 


-32 


-19 


22 


33 






— 














-32 


-19 


-33 


-24 




TT+TT 


38±27 


42 


54 


43 


42 


23 


22 


26 


31 








-49 


-70 


-69 


-85 


-32 


-19 


-72 


-36 


AS = 1 


Exp. 


60° 


90° 


60° 


90° 


60° 


90° 


60° 


90° 






- 


10 


12 


9.5 


11 


1.2 


1.1 


5.5 


5.9 






3±37 


10 


12 


9.5 


11 


1.2 


1.1 


5.5 


5.9 




K+TT- 


-9.3±2.9 


-7.5 


-7.7 


-6.1 


-6.4 


-6.1 


-5.5 


-6.7 


-7.0 








19 


23 


18 


22 


7.4 


7.6 


13 


15 


B+ 


K+tt" 


1±12 


-15 


-15 


-11 


-12 


-6.1 


-5.5 


-9.1 


-9.8 




K°7T+ 


1±6 


3.4 


3.9 


2.9 


3.4 


1.2 


1.1 


1.3 


1.3 




K+K- 




-12 


-16 


-13 


-11 


-6.1 


-5.5 


-6.7 


-7.0 








5.0 


7.9 


4.7 


3.6 


1.2 


1.1 


1.3 


1.3 


Bs 


VsVs 




8.1 


11 


9.5 


8.7 


1.2 


1.1 


3.3 


3.5 














-80 


-71 














TT+TT 




-1.4 


-1.3 


-0.9 


-1.3 












7r°7rO 




-1.4 


-1.3 


-0.9 


-1.3 











Table VI: Fit results for direct CP-asymmetries (Acp). 



Sff- (%) 




SIT(3)- 


claslic 




()ii-slu>ll cc 


I 


II 


III 


IV 


AS = 




Exp. 


60° 


90° 


60° 


90° 


60° 


90° 


60° 90° 




K+K- 






-99 


-83 


-99 


-90 
















-30 


2.3 


-27 


-5.2 


-19 


-10 


-24 -16 




VsVs 






-23 


-18 


-49 


-36 


-19 


-10 


-97 -86 










6.2 


-1.2 


6.2 


-1.2 


-19 


-10 


-24 -16 








-58±34 


-69 


16 


-67 


-1.3 


-81 


-6.5 


-81 -9.7 










-67 


8.4 


-39 


51 


-19 


-10 


65 86 



Table VII: Fit results for the time-dependent CP-asymmetry parameter Sff. 



23 



